function diff = q_0_p_fun(m,param)

    if param.K == 0 
        diff = exp((1-param.alpha).*2./param.sigma^2.*int_delta_OJS(m,param)).*(1-param.alpha).*2./param.sigma^2.*(param.zeta+param.s*param.lambda)./m;    
    else    
        diff = exp((1-param.alpha).*2./param.sigma^2.*int_delta_fun(m,param)).*(1-param.alpha).*2./param.sigma^2.*(param.zeta+param.s*param.lambda)./m;
    end

end